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Abstract. We study a model based on scalar complex fields coupled to a scalar real field, where 
all fields are treated classically as c-numbers. The model describes a composite particle made up 
of N constituents with bare mass niQ interacting both with each other and with themselves via the 
exchange of a particle of mass /Uq- The stationary states of the composite particle are described 
by relativistic Hartree's equations. Since the self-interaction is included, the case of an elementary 
particle is a nontrivial special case of this model. Using an integral transform method we derive 
the exact ground state solution and prove its local stability. The mass of the composite particle is 
calculated as the total energy in the rest frame. For the case of a massless exchange particle the 
mass formula is given in closed form. The mass, as a function of the coupling constant, possesses a 
well pronounced minimum for each value of ^o/mQ, while the absolute minimum occurs at /iq = 0. 



1 Introduction and results 

In this paper we derive and study the exact ground state solution in 3 + 1 dimensions for a 
self-interacting system whose dynamics is governed by the relativistic Hartree's equations. 
Our motivation is to study bound states formed solely by self-interaction. Such bound states 
(sometimes called nontopological solitons) offer a possibility of understanding the internal 
properties of particles, such as their masses, charges and magnetic moments. Self-interaction 
is a purely nonlinear and nonperturbative phenomenon which is neither well understood 
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nor properly appreciated at present. However, we argue that in the quantum domain it 
is rather important and deserves close study. It may be the crucial missing element in 
our understanding of quantum phenomena. Self-interaction appears frequently in quantum 
field theories such as, for example, in quantum electrodynamics, where it leads to infinities. 
Intuitively, the reason why these infinities appear is that we are using a perturbation theory 
in which particles are associated with free fields such as plane waves or wave packets. The 
fact that one is able to cure the theory and remove the infinities in a self-consistent manner 
can be seen as an indication that the original theory, before the perturbation theory based 
on free fields is applied to it, is correct and should merely be treated differently. What is 
needed, perhaps, is an approach which is based on the self-interaction bound states in place 
of free fields. An attempt of such an approach for the case of quantum electrodynamics 
will be presented in a forthcoming paper. In the present paper we attempt to present the 
self- interaction bound states in their own right, within the context of a classical field theory. 
In doing so we will concentrate on their internal properties, in which they differ so much 
from free fields. 

The following model is motivated by its simplicity, self-consistency and by the fact that it 
is connected to the large Nc limit of QCD. It is not intended to be realistic in first place. Yet 
some of its elements, for instance relativistic invariance, dimension 3 + 1 and the nonlinear 
interaction, are clearly realistic. 

We consider a system of N complex scalar fields (r, t), j = 1, . . . , and a real scalar 
field $ (r, t) with the lagrangian 

N 1 

C = -Y: {d'^*d.^j + ml^*^, - g<!>^ {^*^,y) - - (d''<^d.<^ + /ig*^) , (1) 

where p = q = 1, and the equations of motion 

(□-m2 + (7$(r,t)) ^,(r,t) = 0, (2) 

N 

(□ - 111) $ (r, t) = -gY: ^* (r, t) (r, t) , (3) 

i=l 

where □ = A — d"^ /dt^ = d^dy. All the fields are assumed to be classical, i.e. commuting. 
Notice that the interaction term in (|ID is not positive definite and hence the stability question 
arises. In this paper only the local stability is investigated in full detail. Here it suffices to 
say that stable one-particle and many-particle bound states exist, with energies below the 
threshold to the continuum states, provided certain restrictions are placed on the couphng 
constant g and on ^q/ttiq. The fields are normalized according to 

= 1, j = l,...,Ar. (4) 

Notice that the left-hand side of is a constant of motion. In the present context the 
1 on the right-hand side of is purely conventional, in fact we could equally well chose 
any other number. However, equations (0) and allow for a rescaling of the fields, which 
can be chosen to restore the 1 in @. The fact, that are complex fields and hence are 



d^i 



9^* 
dt ^ 



' dt 
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normalized according to (^), is an essential one since otherwise the model would contain the 
case fiQ = mo, \1/ = $ as a special case and then would be identical to the standard $^ model 
which is unstable even locally. 

Equations (H-^ admit plane wave solutions with continuous energies \E]^\ = k^+mg > mg 



vE',(r,t) 



ik, -r—iEi^ t 



$(r,t) 
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N 



n 



(5) 



where Q is an arbitrarily large normalization volume (a detailed discussion of plane wave 
solutions in the context of nontopological solitons is given in ref. As will be shown 

below, the characteristic feature of these states is that for any N their total energy is always 
positive and larger than the corresponding bare mass Nmo, i.e. -Etotai > Nmo and hence on 
the energy scale these states fill the continuum above Nitlq. 

In contrast to the above continuum states there exist bound states with total energies 
< E'totai < NrriQ. In the rest frame of the bound states the corresponding stationary state 
solutions of @ and (|^) are of the form (r, t) = ipj (r) e"*^^* and $ (r, t) = (r) /g, where 
\Ej\ < tjiq. For these stationary fields the equations of motion (j^), 
condition (||) become respectively 



and the normalization 



(A-7j + 0(r)) 



(6) 



(A - (r) 



N 



(7) 



and 



d^i 



2\Ei 



J 



N. 
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where 7^ = — Ej > 0. Once a solution is obtained, it can be Lorentz-boosted to 
an arbitrary inertial frame. Equations and (|^ must be supplemented with suitable 
boundary conditions (they are discussed below). 

This system constitutes a model of a composite particle made up of spinless con- 
stituents of bare mass mo and the constituent energies \Ej\, which interact with each other 
and with themselves via the exchange of a particle of mass fiQ. We identify the physical 
mass M of the composite particle with its total energy -Etotai- The total energy follows in 
the usual way from the lagrangian (|ip 



M 



r ( ^ 

= Etotai = J lj2 [v^; ■ v^, + [e] + ml) 



+ 



29' 



The gradient term (V*/))^ can be eliminated using (|^) and we obtain 



/■ ^ r 1 

M= dhj2 v^; ■ vij, + [e] + ml) r,^^ - -(pr^i^, 
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Eliminating the remaining gradient term by the aid of (j^) and using the normahzation 
condition (||), we get 



M = f:(|E,| + i/d=^r0(r) |^,(r)|^) . (9) 

For the continuum states given in (^) the last term in (P) vanishes as Q ^ oo and hence 
Etotai = J2jLi \Ej\ = YljLi l-^kjl > A^mo- In the case of bound states the last term in (^) 
does not vanish and hence the total mass is not identical to the sum of the constituent 
energies \Ej\ (i.e., to the sum of the bare masses mo and the binding energies mo — \Ej\ of 
the constituents) but contains an additional term. Below we will show that for stationary 
states this term is strictly positive. While the binding energies can be interpreted as the 
energies which are gained by putting particles in an attractive potential, the last term in (^) 
represents the energy which must be spent to create the potential itself. In a consistent field 
theory these two sorts of energies are interdependent and always appear together. 

In the case of = 1 the mass formula (P) defines the mass of an elementary particle, 
which we denote by m. It will be shown below that for stable bound state solutions of (0-0) 
we have m < mo and hence the ground states of particles are bound states and not the 
continuum states for which m = mo. It is therefore not an assumption, but a necessity, 
that in this model the particles are described by self-interaction bound states rather than 
by plane waves or wave packets. The main aim of this paper is to solve (|§-0) for the ground 
state and to determine the dependence of the total mass M and the size parameter tq on 
the coupling constant g, the bare mass mo, the mass of the exchange particle fiQ and the 
number of constituents A^. 

In a situation in which all of the constituent particles are in the same spherically sym- 
metric state with \Ej\ = \E\ < mo and ijjj (r) = (r) = ^p (r) /^/Att, the system of equa- 
tions (H-0) simplifies to 

(^ + '|:-7^ + <^w)^W = o (10) 

and 

^2 d A 9'N , 

( d;:^ + -d;:-/^oj0(r) = -^l^(r)| , (11) 

where 7^ = mg — E"^. Since we are mainly interested in the properties of the ground state, 
which is the lowest energy solution of ([IO| - p]TD , we can restrict ourselves to this situation. For 
the mass we obtain 

M = N (^\E\ + ^J dh(f){r) \ilj{r)f^ . (12) 

Contrary to the case of plane waves, the particles in the present model are extended. As 
a measure of their spatial extension, we define the size parameter 

ro = 2\E\ J d^r r |^(r)|^ , (13) 

where the factor 2\E\ stems from the normalization condition (Bl). 
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The stationary states of the composite particle are described by Hartree's equations. This 
follows from the fact that (0) possesses the well known solution 



2 N 



(r) = 2mV (r) = — ^ y d^r 

i=l 



g-Aio|r-r'| 



/m2 



^.(r')l 



(14) 



which, if inserted in (^), gives the usual Hartree's equations with the effective potential V (r). 
Notice that the self-interaction {i = j) is included in the Hartree's equations and therefore 
the case of an elementary particle, where = 1, is a nontrivial special case of this model. 
Also notice that in the present model the Hartree's equations are exact and not approximate 
as in the usual case. The model would be more realistic if the constituent particles were 
fermions. Yet, despite this shortcoming, it is useful. For instance, it was shown in ref. that 
baryons in the large Nc limit of QCD are described by nonrelativistic Hartree's equations 



+ 2me + (f (r) 



X 







where 



2mg 
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3 / iXir 



d'r' 



(15) 



(16) 



e is the binding energy, m is the quark mass and x is the nonrelativistic single-quark wave 



function which is normalized according to / d r|x(r) 



1. Very briefly the reason for 



this is the following. The color part of the baryonic wave function is a singlet, which is 
antisymmetric, and hence the rest of the wave function must be symmetric. The color 
and spin degrees of freedom can be neglected in the lowest order of approximation and, 
at the scale of hadronic bound states, the linear (confining) part of the potential can also 
be neglected. Neglecting the relativistic kinematics and certain many-particle effects, the 



resulting equations are the nonrelativistic Hartree's equations ([T5| - [16D for identical quarks. 
In this paper we will include the relativistic kinematics and solve the relativistic Hartree's 
equations (^0), since our method works equally well in this case. 

Thus, applied to the quark model, V (r) in (|T^) with /xq = can be interpreted as an effec- 
tive potential in which dressed constituent quarks are bound. The quark-quark potential for 



point-like quarks is g j 



An (2m) 



, which is obtained by substituting (r')| = 6 (r') /2m 
in (0). Comparing this potential with the QCD motivated quark-quark potential for point- 
like quarks, 4as/3r, we obtain 



9^ 



4«s (m) 



(17) 



47r4m2 3 

where ag (m) is the strong running coupling constant of QCD taken at the mass m of 
the constituent quark. Equation (|T^) can be used to relate the results obtained for the 
present model to the large Nc limit of QCD. The relation to the nonrelativistic Hartree's 
equations (|l5|-0) is determined by g''^ = g'^N/ATcAm \E\ and 2me = —7^. Witten drew two 
main conclusions in arriving at Hartree's equations (|15|-|16|). First, provided that g' does 
not depend on Nc (i.e. g oc l/v^iVc), the baryon masses increase linearly with Nc. This is 
an immediate consequence of (0). Second, under the same provision, the size and shape 
of the baryons do not depend on Nc. This follows from (|13D for the size and from ( p!0| - |TT|) 
for the shape. However, to derive further conclusions concerning, for instance, the mass 
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spectrum or the dependence of the mass on the couphng constant, requires the solution of 
the system ([To| - |llD . We will take up these issues here in a more general context where fiQ 
is arbitrary. Besides the connection to QCD, the model is useful on its own right since it 
is self-consistent, the dimension is natural (3 + 1) and the interaction is realistic (exchange 
of particles) and includes the self-interaction. Therefore, it allows us to study the internal 
dynamical properties of particles such as the relation between the bare mass mo and the 
physical mass m, or the dependence of m on the coupling constant g. 
The main results of this paper are as follows. 

Mass. For /iq = there is a doublet of spherically symmetric solutions in the = 1 
sector, which corresponds to a doublet of elementary particles with masses m and m*, where 
m < m*. Consequently the composite particle with constituents possesses N + 1 states 
with spherical symmetry, which can be classified according to how many of the constituents 
are in the excited state. The lowest state of the composite particle on the energy scale is the 
ground state with the mass given by 



N 



Nmt 



/ 














\ATTmlaoJ 







+ 



1 - 



J 



where oq is a numerical constant. 

For Ho ^ there is a triplet of spherically symmetric solutions in the = 1 sector, 
which corresponds to a triplet of elementary particles with masses m, m* and m**, where 
m < m* and m < m** . Consequently the composite particle with constituents possesses 
now A^ (A^ + l)/2 + A^+l states with spherical symmetry, which again can be classified 
according to how many of the constituents occupy one or the other excited states. The 
dependence of M/Nttiq on g'^N/Airml for various fixed values of fio/'^o is illustrated in 
figure |l]a. In this figure the ground state is plotted with a solid line, the excited state, 
having all of the constituents residing in the state *, is plotted with a dashed line and the 
excited state, having all of the constituents residing in the state **, is plotted with a dotted- 
dashed line. Notice that, for each value of /lo/mo, the mass M of the ground state acquires 
a local minimum at the maximally allowed value of g'^/AirmQ (dotted vertical line), which 
becomes the absolute minimum M/Nmo = 2^2/3 in the case /io = (dotted horizontal 
line). The appearance of such well pronounced minima of the mass may be a phenomenon 
which is more general than one restricted to the present model (see ref . , for instance) and 
hence may have some deeper significance (some speculations were discussed in ref. [|]). 

Size. For /Xq = the size parameter of the ground state is given by 



To 



\/25o 4:7iml 



mo 



/ g'N Y 



where Sq and cto are certain numerical constants. For the ground state ro — > oo as g/mo — >• 0. 
This singularity disappears when /zo 7^ 0. For the case of minimal mass the size parameter 
of the ground state becomes minimized too and is given by tq = a/25o/«o"^o- 

Stability. The stability properties are investigated according to three independent cri- 
teria. The first one is based on the mass defect AM = Nm — M , where M is the mass of 
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the composite particle and m denotes the mass of its free constituents. We found that the 
ground state is stable against disintegration, i.e. AM > 0. 

If the mass of a particle is larger than its bare mass, then such a particle is unstable. The 
reason is that, in this case, some of its continuum states (for instance the states in (^) are 
energetically preferable to the bound state of the particle. According to this second criterion, 
in the case /zq = the ground state is stable. For /iq 7^ the ground state is unstable for 
smaller values of the coupling constant and stable for larger values (see figure |I|a). For 
yUo > 0.36mo the ground state becomes unstable for all possible values of the coupling 
constant. 

The third stability criterion is based on S'^Etotai > 0, which means that a bound state 
solution is locally stable if it corresponds to a local minimum of the total energy {6Eiota\ = 0). 
This criterion was shown by Rosen 0] to be necessary and sufficient for a dynamical stability 
in the sense of Liapunov. In the present paper we prove that for /ig < y^uiQ the ground 
state is locally stable. Moreover, it is proved that our model is the only option with locally 
stable bound states among the class of theories based on lagrangians of the form (|l|), where 
p and q are arbitrary real positive constants. 

Existence condition. A short range interaction cannot produce a bound state unless 
the strength of the interaction is sufficiently large, i.e. the coupling constant is larger than 
a certain limit. If the system is relativistic then, in addition to this constraint, the coupling 
constant must be smaller than a certain limit. For the present model this existence condition 
is given by 

N AT[ml ~ Auml ~ N Aiiml ' 

where go and gi are certain functions of fxo/^o and are independent of A^. 

Wave functions. Equations (|T0|) and (pT]) are solved using an integral transform method. 
The resulting wave functions are plotted in figure ^ and for various values of the 
parameter a = /io/27, where 7^ = mg — E^ > 0. More precisely, the scaled functions 
1^1 1/2 ^ (r) /7'^/^and (r) /7^ are plotted, since in this form the wave functions are both di- 
mensionless and depend on the dimensionless variables jr and a only. Notice that (r) > 
for all r. There are no bound state solutions for the case (r) < 0. 



2 Integral transform 

In ref. we discussed a method for solving the Schrodinger equation, which is most suitable 
if at large distances the potential is proportional to an exponential function. It follows 
from (^) that this is the case in the present situation and hence, according to ref. the 
integral transforms appropriate for (|ToD and ( ITll) are 

POO 

^ (r) = ao e-^'' / d/i e'"'' g U) (18) 
Jo 

and 

0(r) = boe-''°' / d/x e-'^'-y (/i) , (19) 
Jo 
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where is a normalization constant and 60 is a constant. Substituting (|T^) into (|TOD and 
(ITqI) into firri), we obtain 



d/i' 



Mo 



/i' (/x' + 27) d/i' 



y / djj," V (yu' - /i") g (/i" - /io) 

"'Mo 



(20) 



and 



V(/x) = l-^^^(/x + /xo-27) 



d/i' 
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27-Mo /i' (/i' + 2/io) d/j,' 



X 



27-Mo 



d/iX/i'-/i")f?(/i" + /xo-27), (21) 



where 6 [fi — a) is a step-function which vanishes for fi < a, equals 1/2 for fi = a and is 
equal to 1 otherwise and a = /to or a = 27 — /io respectively. In deriving (^) we have 
assumed that /io/27 < 1. It will be shown below that this implies that we are restricting 
ourselves to the case /io < 2mo. Now we introduce the dimensionless variables 



/i /io bo 

'=2^' "=2:^' ^=2^ 



c 



and redefine the functions g and V in terms of these variables 

Ra iv, C, s) = Q (/^) , V„ (r/, C, s) = V (/i) . 
Equations (^) and (|2lD become 

dt d /•* 



(22) 



(23) 



{r^X,s) = l-r^e{s- a) j^^^ - dt' V. (r/, C, t - t') (r/, C, t' - a) (24) 



and 



V.(r/,C,s) = l--^(s + ^-l) / 



dt 



d 



t {t + 2a) dt 



f dt'R,{r]X,t-t')R^{r],C,t' + a-l) (25) 



respectively. Substituting (pSj) in (|^) and performing an integration by parts, the depen- 
dence on Va {rj, (, s) is eliminated and we obtain 



RAv,C,s) = i-vHs-(r) I —-^RAvXt-cx)+ce{s-i) 



dt 



X / dt 



t{t 

,RAvX,t-t') d 



t{t + l) 



t'2 - (t2 



dt' 



dt" R,{7],C,t'-t") R^{r]X,t"-l) . (26) 



We return now to the question of boundary conditions for (|^) and (0) or, equivalently, for ([T0| ) 
and (|TT1) . Since both ( p!0D and (|TT|) are second order equations, they must be supplemented 
by two conditions each. We first consider the case of ip (r) where, as usual, the conditions 
are: ip (00) = and {ip (0) | < 00. The first condition has been taken into account by (plsl) . 



8 



As for the second one, in order to translate it to a suitable and practical condition, we 
write ([T8|) using (^) and (|23|) as 



roo 

^ (r) = 27ao e"^'' / ds e-'^^' R„ (r/, C, s) 
Jo 

and replace {ip (0) | < oo by 



Jo 



< A 



(27) 



(28) 



where A is some suitable real number. This condition becomes the true boundary condition 
as So ^ C)0 and then tq — 0, but for a finite accuracy result some finite values of Sq and Vq 
are sufficient. The order of the limits cannot be interchanged since R^j [r], (, s) as a function 
of s does not vanish at infinity (its behaviour in the large s regime can be described roughly 
as s"cos(6s + c) with some constants a, b and c). The conditions for (r) are the same: 
(oo) = and |0 (0) | < oo. Again, the first condition has been taken into account by (|T^, 
whereas for the second it follows from 



POO 

Jo 



that 10 (0) I < 00 can be replaced by 



ds e-2^^»^ V. (r/, C, s; 



< A 



(29) 



(30) 



Equa- 



with some finite values for tq and sq depending on the accuracy to be achieved, 
tions ( PBD and (^) show that the constants rj, ( and a are not independent. 

Thus, our next task is to solve ( ^6]) with parameters r], ( and a satisfying (pHf ) and (|30D. 
This will be achieved in two steps: first we solve (|26|) regarding rj, ( and a as independent 
and then for each given a we determine rj and ( satisfying (^) and (pO|). Thus, at this 
second step rj and ( become functions of a. 



3 Solution of integral equation 

To obtain the solution of ( ^6|) we make the Ansatz 

[s/(t] [s—na] 

n=0 m=0 

00 00 

= ^ ^ {-ilf C^e{s -na - m) Lpnm {s - na - m, a) , (31) 

n=0 m=0 

where a number in square brackets represents the largest integer which is smaller than or 
equal to that number. Notice that, for a fixed s < 00, the right-hand side of the above 
equation is a finite sum, so the question of convergence does not appear. Substituting ( PT| ) 
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in ( p6D and comparing equal powers of t] and (, we obtain the recurrence relation for the 
functions (Pnm 



° (t + ncr + m + i 



2/ 4 



+ E E / dt f;-^:";'^./'" E E 1; / dtV.-M-. (t' - t'\ a) (t", , (32) 
k=o 1=1 -^0 [f + ka + l) -a^ ,=i / 



i=0 j = 

where n = 0, 1, . . . , m = 0, 1, . . . and 

V^-i.m (y, 0-) = , (^00 (y, 0-) = 1 , (33) 

which guarantees that the Ansatz is in fact a solution. Here and in what follows we adopted 
the convention that if the upper bound of a sum is less than the lower bound then the 
contribution of the sum is zero. Notice that for fixed n and m the evaluation of the right- 
hand side does not require knowledge of the left-hand side. Also notice that the Ansatz 
has been chosen so that the recurrence relation does not depend on the parameters r] and 
C,. Therefore, once the ipnm functions are calculated, (^) is solved for any rj and C, and we 
can use this solution to find the values of 7] and C satisfying the boundary conditions (|28| ) 
and (|30|) . Equation ( pl|) combined with the above recurrence relation constitutes the exact 
analytical solution of (P§). 



The recurrence relation (|3^ can be iterated, so that finally no ipnm functions will appear 



on the right-hand side. For instance, 

(g + l){y + a) 1 [V dt (t + 1 - a) (1 + a) 

^^°^^'"^=^" a{y + a + l) ' = 2^ io (t + 1) (t + 2) (t + 1 + a) (1 - a) ' 

(34) 

However, for higher n and m one obtains multiple-integral representations of ipnm which 
are not useful for our purposes. Since one cannot express the functions ipnm in terms of 
elementary or special functions, one has to devise a numerical algorithm for their evaluation. 
A suitable algorithm is given in the appendix. As a by-product of this algorithm, we also 
obtain a function fnm (s, a), which is very useful since Va (rj, (, s) is related to it in a simple 
way: 

■1 oo oo 

VAvX,s) = l--Y.T. i-vTreis + a-na-m) F^m {s + a - na - m, a) , (35) 

^ n=0 m=l 

where 

Fnm [y, Cr) = / dt fnm [t, a) . (36) 

Jo 

Notice that, for a fixed s < oo, the right-hand side of (^) is a finite sum. 

Once the functions (pnm {y, o") and Fnm {y, c") are calculated we can use them to calculate 
the function R„ (rj, (, s) according to (^) and the function {rj, (, s) according to (^), and 
to determine the values of r] and ( for which the boundary conditions (^) and (^) are 
fulfilled. A typical result is illustrated in figure H. 
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4 Wave functions and binding energies 



To determine the wave function ip (r) of a constituent we have to determine the normahzation 
constant oq. Using the normahzation condition (|^) and equation (p7D, we obtain 



27C 
a\E\ 



(37) 



where 



«=16C/ dee'^'(0, ^.(0 = e-^/ dse-2«^i?.(r/,C,s) . (38) 

Notice that a is the only parameter upon which a depends. Using (pSj) and the function 
Ra iVy C) -s) calculated above, a can be determined. A useful fit is 



a = + aia + a2(J , cr < 1 



(39) 



where ao = 3.52(2), ai = 10.9(2) and ^2 = 3.82(9). Using (jS^) and (|7D we can write the 
wave function in a dimensionless form: 



V a "'0 



7' 



3/2 



(40) 



Notice that the right-hand side of (^) is a function of the dimensionless quantities jr and 
0" only. Similarly, using ( ^9]) and (|22D the potential function (r) also can be written in the 
dimensionless form 



1 



47] e-^"^' ds e-2-^" V. (?7, C, s) 

JO 



(41) 



In figures |^a and ^ we plotted the wave function and the potential function according to 
and (^If ) respectively for various values of a. 

We define the binding energy E — rriQ of a constituent particle as the difference between 
its constituent energy E and its bare mass tuq. To determine the constituent energy E we 
substitute (|37|) into the expression for ( provided by (|2^) and using 7^ = ml — E"^ obtain 



EUlml-E 



1 g^N 

2 47ra 



Solving for \E\ we obtain two solutions: 



(42) 



mo 



1 



\ 



'_IN_ 



and 



mo 



1 + 



(43) 



(44) 
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which are real if g'^N/Airmla < 1. The first root is the constituent energy of each constituent 
in the excited state, while the second root corresponds partly to the ground state and partly 
to the excited state, depending on the value of a. Equations (E3|) and imply that 



mo V ^ 

and 

yZ mo 

Notice that e* is not larger than e, and nevertheless e* is assigned entirely to the excited 
state. The reason is that, as will be shown below, the total energy (the mass) corresponding 
to the bound state with constituent energy e* is always larger than in the case of a bound 
state with constituent energy e. 

In the case when /xq = 0, the function a = ao = 3.52(2) is a pure numerical constant 
and hence (|4^) and (|^) are explicit formulae. Otherwise, a is a function of cr which itself 
is a function of both parameters fio/mo and g. We will not need the explicit form of this 
function, since a will serve only as a parameterization variable at fixed ^q/uiq. The allowed 
range of a follows from 



a = = , 47 



which is obtained from the definition a = /io/27 and from 7^ = rriQ — E"^ . Substituting (|i5| ) 
and (^61) in (0), we obtain respectively 

< , < 1 i^i (48) 

2 mo V2 mo 

and 

45£^<-<oo. (49) 
The coupling constant g'^N/Airml, for instance, is parameterized by a at fixed ^o/mQ via 



pL = ^Jl-^-(J^)\ (50) 

Equation ( |50D is obtained if one combines ( ^2]) and ( ^Tj) to eliminate i?. Using (|50| ) we can 
convert the dependence on a in ( ^3|) and (|^) to a dependence on g^N/inml and /xo/mo. 
The result, e/mo and e*/mo as functions of g'^N/ATiml and fio/rrio, is plotted in figure^. The 
dotted line indicates the value l/\/2 which is the boundary value deviding the two functions. 



5 Existence condition 

Equations (^) and (Q) imply that 



< 



a 

< — 

- N 



(51) 
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since otherwise the constituent energy becomes a complex number, in which case the bound 
state does not exist. In the case when /zq = 0, a = = 3.52(2) is a pure numerical constant 
and hence (|5TD is an explicit necessary and sufficient condition for the existence of bound 
states. 

In the case when hq ^ (|5T| ) is not suitable as an explicit existence condition, since 
in this case a is not a pure number but rather a function of a and thus of g'^N/AirmQ and 
fio/mQ. To obtain the existence condition for the case when /xq 7^ 0, we write (0) as 



9' 




(52) 



which is the parametrization of the coupling constant in terms of a at fixed /io/Trio and N. 
The inequality on the right-hand side of (|5^ ) becomes an equality for a = yUo/v^^o- Since 
a varies in a definite range given in (^) and (^), equation ( p2D determines the allowed 
range of the coupling constant, in which a bound state solution is possible. In figure ^ we 
plot g'^N/AiTml as a function of a for a few typical values of fio/rrio (solid lines). Notice that 
g^N/4:7Tml possesses a local maximum in the region of smaller a and a local minimum in 
the region of larger a. From (|5^ ) and ( p9D we obtain 



c"i = ^ I 1 ; 

V2mo 8ao \moJ 

for the location cr = o"i of the local maximum and 



l^^^(l!2.)\ol(ilY] (53) 




.^,/-.0((^)^l (54) 




for the location a = of the local minimum. The respective accuracies of (^) and (Q) are 
sufficient for all practical purposes. This can be judged by looking at figure |^, where ui from 
( [53| ) is plotted as a vertical dotted line and do = \J 010/0.2 is plotted as a vertical dashed line. 
The maximum and the minimum split the allowed range of g'^/AnmQ into three regions: one 
to the left of the maximal value 



9^ < 1 ^1 li^ 
AirmQ ~ N Airml ' 2 mo 

one between the maximal value and the minimal value 



< — — 2 < ~7Z — 2 ' o — < o" < o"! 5 (55) 



^0/^^/1 9f 



< 2 - 2' '^i<o"<o"o, (56) 



and one to the right of the minimal value 



^-^<-^, ao<a, (57) 



where 



47rmo o"imo V 4 Vaimo/ v2mo VVmo/ / 
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= a (a,) J^Jl - i (J^y = (2^+ aO ^ + O ((1^)") (59) 



and a is given in (pQl) . These three branches of the couphng constant correspond to three 
different states (particles), and below it will be shown that (p6D corresponds to the ground 
state while (^) and (|57D each correspond to a different excited state. Hence for the case of 
/io 7^ the model predicts a triplet of particles. Notice that in the case when /xq = ( |5BD 
becomes 

0<7A<^- (60) 

If for certain fio/rriQ and it happens that g'^/Anml is outside of the bounds defined in (|56|) 
then a ground state solution and a particle associated with it do not exist and there are no 
stable particles in that case. Thus (|56D is an explicit condition for the existence of a ground 
state and its associated particle. Notice that (|56| ) is in qualitative agreement with the 
nonrelativistic case of the Yukawa potential , where the ground state solution is possible 
only if the coupling constant /^/47r is larger than (1.6798/2) (/xo/mo). 

In the case when the condition (|56D for the existence of the ground state is violated 
while ( |55D is fulfilled, the particle associated with the ground state does not exist while the 
excited state does. However, as will be shown below, the mass of such a state is larger than 
its bare mass. Hence this state is unstable and represents an unstable particle. 



6 Mass and bare mass 



Instead of calculating the last term in (|T2|) directly, we can use the virial theorem (this 
theorem is proved in section ^ below) part of which is the relation 



d^r (f){r)\ip (r)| 



2 mo /mo \E\ 



+ 



1 fil 



\E\ mo J 3g^N 



J dh 02 (r) 



(61) 



Substituting ([61| ) in (|T2]) we obtain 



.r /2|E| Imo 1 /i^ 



3 mo 3 I -El SniQg^N 



dh (J)'^(r)] > . 



(62) 



Notice that the total energy ( p^D is positive definite. 

We first consider the case when fiQ = 0. In this case the last term in (^2[) vanishes and 
we obtain 



2 /lEl 1 mo \ 



Substitution of (|1) and 
V2 



m 



3 V^o 2 \E\ 
yields 



(63) 



/ 



N 



Nmn 



N V 



[^TrmlaoJ ^ g^N \^ y 1^47rm2aoj 



(64) 
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for the ground state and 



M, 



TV AT 



Nmr 



V2 



1 




[i 









1 - 



1 + 



(65) 



for the excited state. The masses of the corresponding constituent particles are 

Ml, m* = Mil ■ 



m 



(66) 



Thus, for the case of /xq 
(A^ = 1) with masses m, 
and Ntuq 



the model predicts a doublet in the elementary particle sector 
I*. Notice that (p^|66D imply that for fixed parameters g'^N/A 



Mn < M, 



NN 



m < m* 



TTTTln 



(67) 



despite the fact that according to (^) and ( ^31) e* < e. The reason is that, besides the 
contribution from the constituent energies (first term in (^), the mass receives also a con- 
tribution from the energy (second term in (^) which must be spent to create the potential 
in which the constituents are bound. This contribution compensates the difference between 
e* and e and places M^n above M^- Also notice that while Mnn oo and m* ^ oo as 
g 0, Mn and m stay finite. 

The state with mass Mnn is the highest excited state with spherical symmetry, since 
we have assumed that all of the constituents occupy the same state and that there are only 
two bound states which each of the constituents can occupy. Our method, however, can be 
extended easily to a situation where one or more constituents occupy the excited state with 
the corresponding constituent energy e* while others remain in the other bound state. The 
resulting excited states will possess masses between Mn and Mnn- 



Mn < M^nM: 



IN, 1VJ-2N, 



< M 



NN 



(68) 



These are all the possible spherically symmetric states, since otherwise there would be more 
states than just a doublet in the A^ = 1 sector. Hence for the case of /xq = 0, the total 
number of states with spherical symmetry is 



A^+ 1. 



(69) 



A further consequence of (|6^|66D is 

AMtv = Nm - Mn>0, 

which means that the ground state of the composite particle is stable against disintegration. 
Similarly, one obtains 

AMnn = Nm* - Mnn > , (70) 

which means that the excited state cannot decay by disintegration with emission of particles 
of mass m*. 

Above we discussed the stability of the composite particle from the point of view of 
disintegration. We now discuss the stability from the point of view of a comparison between 
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the bare mass and the physical mass. This criterion allows us to discuss the stability of even 

I it follows that 



the elementary particles of the theory. From 



< 



9' 



< 



N 



Mn < NrriQ , 



(71) 



for the ground state and 



< 



< 



x/3 



2 N 



< 



9 



-\/3 c^o 



" N 



^ Mnn > NniQ 



^ Mnn < Nmo 



(72) 



(73) 



for the excited state. If we want to associate a particle with a self-interaction bound state, 
then the result that the physical mass is larger than the bare mass means that such a particle 
is unstable. The reason is that in this case some of its continuum states are energetically 
preferable to the bound state. Thus from ( ffl]) we conclude that the ground state is 
always stable if it exists, while the excited state Mnn is unstable for smaller values of the 
coupling constant g in the case of ([7^), and stable for larger values of the coupling constant 
in the case of ([731). The phenomenon of instability due to the mass being larger than the 
bare mass is not new. It was observed, for example, for the ground state in the model of 
ref. 0, which is based on the nonlinear Schrodinger equation in 3-D. In that model one can 
give a formal proof (see ref. for the details) that the ground state is unstable in the sense 
that its wave function suffers a collapse; i.e., being initially smooth it develops a singularity 
within a finite period of time when subjected to a small perturbation. This means that such 
a state cannot be associated with a stable free particle. 

In the case of ([73| ) the composite particle can exist as a stable (with respect to the 
continuum states) particle with the mass M^v or M^vat. For the mass difference we have 



< Mnn -Mn < Mn 



(74) 



which amounts to less than 6 %. 

We now discuss the case when /iq ^ 0. Using (p!9D , (^2)), (^2)) and (|47|), we obtain 



2 m^g'^N 



d^r 0^ (r) 



'mo _ \E\ \ 
mo) 



(75) 



where 



16(T^?7^ 



a 



roo roc 

Jo Jo 



(76) 



and Vcr (?7, C, s) is defined in (pSf). Notice that /5 is a function of a only. To get some idea of 
the behavior of /3 the following fit is useful 



a < 1 



(77) 
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where /3o = 0.834, /5i = 0.972 and /?2 = 2.93. The accuracy of (^) is about 1%. Using (|7|) 
we can rewrite (^2|) as 



M = Nmo: 



mo 2 \L\ 



(78) 



Substituting (|43|) and (|^) in ( ffSD we can obtain similar mass formulae to (164^661) . More 
convenient formulae are obtained if, using (]5^), we eliminate g'^N/Anml from ( ^3D and (^) 



and substitute the result in ([Tq) . The resulting mass formulae are 



Nrrin 



mo 



V2mo 8ao Vmo/ 




Nrrin 



mo 



(Ti > cr > 



1 /ip 

2 mo 



(79) 
(80) 



and 



A^mo S ( cr, I , oo > cr > ctq 
mo/ 




« (-)■ 



where 



cr. 



mo 



6 \ crmo 
4 V<^™'0, 



'i - i f^V 

4 \amoJ 



1- P / flQ 



6 



\ yUo 47rmQa 



(82) 



The advantage of (|79| - Fll) is that, for a fixed /io/mo and NniQ, they depend solely on a. The 
dependence of M^r, Mq^n and MmN on cr for various fixed values of fio/niQ is illustrated in 
figure 03. In this figure Mjsi /Nmo is plotted with the solid line, Mqj^n /Nmo with the dashed 
line and M^on /Nmo with the dotted-dashed line. The dependence on a can be converted 
to the dependence on g^N/inml via (|52D 



a 



1^0 

amo 



4 Vcrmo 



(83) 



The result is illustrated in figure |l|a. In both figures the dotted vertical lines correspond to 
the maximal values of g'^N/AT^m^ according to ( [52| ) (cr = cii), while the dotted horizontal 
line corresponds to the absolute minimum of the mass M/Nmo = 2^2/3, which occurs 
at /io = 0. Hence we observe that the mass Mjv of the ground state acquires a minimum 
value at the maximal value of g'^/AirmQ. We make a similar observation with respect to 
the maximal value of the mass M^v: this occurs at the minimal value of g'^/ATcm'^ according 
to ( |52D (cr = cTo). Notice that, except for the case when fio/^o = 0, part of including 
the maximum resides within the continuum, i.e. above Nmo, signaling an unstable state. 
This part becomes larger with increasing fio/mo and starting with fio/mo = 0.36 all of the 
mass Mn resides within the continuum. This means that for a stable ground state the mass 
of the exchange particle cannot be too large: /io < 0.36mo. 
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In the case of elementary particles of the model {N = 1) we have a triplet with the 
corresponding masses 



m = Ml , m* = Moil , m** = M^i , (84) 



which are obtained by putting = 1 in ([79| - |8lD . This fact is in contrast to the jj,^ = 



case where there was a doublet, and is a clear indication that the mass is not an analytic 
function of /xq at /iq = 0. The mass is also nonanalytic at g = gi, where it has a minimum, 
and at g = go where it has a local maximum (see equations (^) and ([5P|)). Otherwise it is 
an analytic function of g'^N/AnrriQ and fiQ/rriQ. 

The appearance of a triplet with different masses and different constituent energies at 
the same values of g'^N/ATTrriQ and ^o/mo indicates that there are other states with spherical 
symmetry, which can be labelled by the number of constituents in each of the two excited 
states. The lowest mass state M^r is the ground state, where none of the constituents reside 
in an excited state. Then there are states Mk^n-k,N where n constituents reside in an excited 
state, k of them in one of the excited states and n — A; in the other one. Since all these states 
are spherically symmetric, our method can be easily extended to include them. Basically, 
the only significant change in this case is the tripling of the basic equations, so that we would 
have six coupled equations instead of two. Since in the present case there are no other states 
beyond the triplet, it is fairly clear that these are all spherically symmetric states of the 
model and that their total number is 

A^(A^ + 1) 

iVstates = + N+1. (85) 

Thus in this model there are 3 elementary particles (triplet), 6 composite particles with 2 
constituents (sextet), 10 composite particles with 3 constituents (decuplet) etc. 

One of the most basic properties of the masses as functions of g'^N/Anml at fixed /io/^^o 

is 

i!iI<Z^ ^ JL, (86) 

inniQ ~ 47rmo Nttiq ~ NniQ ' 

where M' and M are any of the masses Mat, M^on, Monn, ^* and m** evaluated at 
g''^N/47rmQ and g'^N/ATrrriQ respectively. An immediate implication of this result is that 

AMn = Nm- Mn > 0, (87) 

AMoNN = Nm* - MoNN > , (88) 

AMnon = Nm** - Mnon > . (89) 

Thus none of the composite particles can decay by disintegration and in particular the ground 
state is absolutely stable whenever its mass is below the bare mass. In the case when any of 
the masses is larger than the corresponding bare mass, for instance M^qn > Nmo always, 
the particle with that mass is unstable as discussed above. 

One of the striking features of the mass (total energy) as a function of the coupling 
constant g'^/Anml is the appearance of a strongly pronounced unique minimum at each 
fixed value of fio/'^o and N (see figure |T]a). These minimal values of the mass and the 
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corresponding coupling constants can be determined from (|79D, (0) and (|53D and, to first 
order in iiq/uiq, are given by 

where ao = 3.52(2), ai = 10.9(2), j3o = 0.834. The fact that the physical mass of a particle 
possesses such an absolute minimum was noticed earlier in ref. [^] in connection with some 
other model, where it was emphasized that this phenomenon might be a general one and 
might have some deep significance (some speculations are discussed in ref. [|]). 



7 Size parameter 



The size parameter ro has been defined in (|T3|). Substituting (^) and (^) in (0), we obtain 

' (91) 



6 

ro= - 



7 ^ml - 



where 



16C 



a Jo 



POO 

Jo 



(92) 



and (0 is defined in (|38D . As in the case of a and P, the function 6 can also be easily 
computed, once the integral equation ( PUD above is solved. A useful fit, which is accurate to 
about 1%, is 



a ^ 



a <1, 

where 6o = 8.4, 6i = 1.35, 62 = 5.75 and a is given in (|39|) . Thus 

1.3 < 5 < 2.4 . 

In the case when /io = 0, we have a = and hence using (^) and (|4^), we obtain 



ro 



S ^_ V260 Anml 
ee* mo g'^N \^ \ 



^47rmgao, 



(93) 
(94) 

(95) 



for the size of the ground state and 



mo g'^N \ ^ \ATTmlao) 



(96) 



for the size of the excited state. Notice that rg ^ 00 as g/rriQ — > while stays finite. This 
singularity disappears when /iq 7^ 0, since in this case g > go > (see equation (pUD). For 
the case of minimal masses, equation ( pO|) for fiQ = 0, the sizes are 

. 1 72^0 



mo ao 



(97) 



For the ground state ( ^7\) corresponds to the minimal size, while for the excited state - to 
the maximal size. 



19 



8 Derrick's theorem and proof of local stability 



Derrick's theorem refers to time-independent solutions of a class of nonlinear equations 
for real scalar fields and it consists of two parts. The first part is the virial theorem and the 
second part is the proof of local instability. We now recall Derrick's theorem and discuss some 
of the problems connected to it. We then generalize this theorem to a class of theories which 
contains the model investigated in this paper, and derive the conditions of local stability. In 
particular we prove the local stability of the ground state solution studied in this paper. 
Consider the lagrangian for a real scalar field 9 



^=o hr -{-^df-fiO)] (9^ 




and the corresponding equation of motion 

A^^-|? = ^/'W, (99) 

where / is a smooth function. A time-independent solution 6 (r) of (^) corresponds to the 
extremum 5H = of the energy functional 



H 



J dh[{V9 f + f{9)] =h + l2. (100) 



Using this fact and a particular form of the variation 5H, Derrick proved that the kinetic 
part Ji and the potential part I2 are related according to 

/i + 3/2 = 0. (101) 

Equation ( |101| ) constitutes the virial theorem. In the case f [9) > this theorem precludes 
the existence of time-independent solutions of (p9|) since in this case both I2 > and /i > 0, 
which contradicts ( p.01| ). If f [9) > is not valid, the energy H is not bounded from below 
and hence a time-independent solution of (^9]) can be stable at most locally. However, using 
a particular form of the variation. Derrick showed that 

6^H = -2Ii<0. (102) 

Local stability requires S^H > (local minimum of the total energy) for all possible varia- 
tions, but to prove the local instability it is sufficient to show that 6^H < for a particular 
variation, so that (|102|) implies that all time-independent solutions of (p9|) are locally unsta- 



ble. Equation ( [L02| ) constitutes the second part of Derrick's theorem. Shortly after Derrick's 
paper Rosen 0] proved that 6^H > is the necessary and sufficient condition for dynamical 
stability in the sense of Liapunov. 

In the subsequent repetitions of Derrick's theorem (see for instance ref. 0) the second 
part of Derrick's theorem was dropped and the local stability condition S^H > replaced by 
the much stronger condition f (9) > 0. However, if one admits time-dependent solutions, in 
particular stationary bound state solutions of the form \1/ (r, t) = e'^^^ip then the latter 
condition cannot be justified in general and the former, contrary to ( |102|) , can be proved now 
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for particular cases. Thus in particular cases at least locally stable bound state solutions are 
possible. A simple illustration of this fact is provided by 



(n-ml + gq{^*^y~^)^ = 0, 1 < g < 5/3 



(103) 



where □ = A — d"^ /dt^. The local stability of stationary bound states of ( |103|) was proved 
in ref. [jlO|- Notice that the energy corresponding to (|103|) is 



H 



d^r 



dm* d^! 



(104) 



and that the interaction potential (last term in (|104| )) is strictly negative {g is a positive 
constant) and unbounded, and nevertheless the ground state has a finite energy and is 



locally stable. The proof of local stability for the nonrelativistic analog of ( |103|) and other 
examples of stable theories with a nonpositive interaction potential can be found in ref. 
and references therein. 

Consider now a system of complex scalar fields {r,t), j = 1, . . . ,N and a real scalar 
field $ (r, t) with the lagrangian 



N 



and the equations of motion 



□ - mo + gq^P 



N 



= 



□ 



$ = -(^p$P-i^(x[/*x[/, 



(105) 



(106) 



(107) 



where mo, g, p and q are real positive constants and the fields are normalized according 
to 



d=^i 



1, J = l, 



N. 



(108) 



dt ' ' dt 

Notice that the left-hand side of ( |108| ) is a constant of motion. There are four important 
special cases to notice: p = q = 1 gives the model investigated in this paper, p = yields 

3 yields the standard field theory, while q 



the case of ( |103| ), and q = 0, p = 3 yields the standard field theory, while q = 0, p = A 
yields the standard field theory. The local stability condition for the second case has been 
derived in ref. [|TU| and the existence of stable bound state solutions in the latter two cases 
has been ruled out already by Derrick's theorem. We now generalize Derrick's theorem to 
the class of theories characterized by ( |105| - |108| ) for the case of time-dependent but stationary 
bound state solutions of the form \E'j (r, t) = ipj (r) e"*'^-'*, $ (r, t) = cf) (r) and for the case 
of arbitrary p and q. Since this class is less general than in the case of Derrick's theorem, 
we will benefit by being able to derive two additional virial relations. 

For stationary fields the equations of motion (|106|) , (|107|) and the normalization condi- 
tion ( |108| ) become respectively 



A 



l] + gq<f{r,^,)''')t 







(109) 
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(A- 



N 



i=l 



and 



1, J 



N, 



(110) 
(111) 



where 7^ 



rrin 



Ej > 0. The corresponding total energy is 



r I ^ 



+ 



(V0)^ + /ig02 



= Hi + H2 — Hj + + H4 
Notice that (|111|) imphes 



Ho 



N 

E 







f mo 








mo /_ 



(112) 



(113) 



In order to make the variation of the energy 6H we have to choose a proper energy 
functional Equations (|111|) tell us that we are dealing with a constrained system. According 
to the standard rules of quantum mechanics there are two ways to perform variation of the 
energy 6H for a constrained system. The first |jll| is to use an energy functional which does 
not depend on the norm of the variational functions corresponding to ipj. The second 
is to introduce the normalization condition by real Lagrange multipliers. Both ways are 
equivalent, but for our purposes the first is more convenient. The unique functional of the 
variational fields ip\j and 0a, which fulfills the above requirement and which reduces to (|112|) 
in the case ipxj = ipj, (j)x = (f) is 



H{X)= Jd'rEU 



-i/d^r [(V0a)' + 



2^2 



(114) 
(115) 



Now to perform the variation 5H we have to choose a set of variational fields ipxj and (px, 
which is large enough to yield all the solutions of 5H = or, equivalently, the equations of 
motion (|109|) and (|110|) . A suitable set of variational fields is 



^IJxj (r) = AV (Ar) , rx, (r) = A*>* (Ar) , <f>x (r) = A^0 (Ar) 



:ii6) 



where t is an arbitrary complex and A and s are arbitrary real numbers. The variational 
fields ipxj are not normalized except for A = 1 in which case the normalization is defined 
in (|111| ). For A = 1 we have 



^1, (r) = V^,- (r) , rij{r) = rj{r) , 0i (r) = (r) , H (1) = H . 



The variations Sipj, Sip*, S(f), 6H and 6^H are defined by 



Sip = 



dX 



tip{r) + r-Vip{r) , 5ip* = 



'dip 



A=l 



dA 



flp* 



v-\7ip* 



:ii7) 



:ii8i 



A=l 
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and 

56 



dA 



s (r) + r ■ V0 (r) , 5H 



'dHjxy 

dA 



6^H 



dA2 



A=l \ / A=l \ / A=l 

(119) 

Notice that at each r the variations Sipj, Sip* and 6(j) are arbitrary and independent of each 
other. Therefore, applying the standard variational procedure of quantum mechanics, we 
conclude that 6H = is equivalent to (|109[ - |110D . 

Since the integration in ( p.l5| ) is over the entire space, we can eliminate the dependence 
of the fields on Ar by making a suitable rescaling of the integration variable and obtain 



H (A) = X'^Hi + H2- X'P+^^'i-^^Hj + X'^'-^Hs + X'^'-^H. 



4 5 



(120) 



where Hi, . . . and Hf, which are independent of A, are defined in ( |112| ). From ( |120| ) and 
6H = it follows that 



2Hi -[sp + 3{q- 1)] Hi + (2s - 1) H3 + {2s - 3) H^ = ^ 
Since ( |121|) must be valid for all s, it implies two separate virial relations 

2H3 + 2Hi -pHj = 

and 

0. 



(121) 

(122) 
(123) 



2Hi -H3-3H4-3{q- 1) Hj 

Moreover, multiplying ( |109|) by ip*, integrating over the entire space and summing over all 
j = 1, . . . , A^, we obtain a third independent virial relation 



N 



Hi + H2-qHj = Y.\Ej 



(124) 



Combining these relations one can obtain other useful relations. For instance, eliminating 
Hi and Hj and using (|113|) we obtain 



N / 

[2 (3 - g) - p] H3 + [2 (3 - g) - 3p] H, = moJ2 i 

j=i \ 



( mo 






mo J 



(125) 



Since H3, H/i and the right-hand side of ( |125| ) are positive numbers, we obtain a necessary 
condition for the existence of bound states 



2(3-g) - p > 0. 
Another useful relation is obtained if we eliminate H^ from 

N 



l[2{3 -q)-p]Hj=moY.^I^ 



mo 



by means of 
+ 2pHi . 



(126) 



(127) 



This is the equation ( |6TD which we used above to determine the mass. 
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Now we compute the second variation of the energy 5^H and derive the local stability 
conditions which ensure that 5^H > for all values of s and t, which are the parameters 
spanning the set of variations. Equation ( |12CI|) implies 



6^H = 2Hi - [sp + 3 (g - 1)] [sp + 3 (g - 1) - 1] if/ + (2s - 1) (2s - 2) H3 

+ (2s - 3) (2s - 4) . (128) 

Eliminating Hi, and Hi by means of the virial relations (|122|) , ( |123|) and (|125|) , we obtain 
^2^ _ 2p (2 - p) s^ - 4p (3g - 2) s + 3p - 6 (g - 1) (3g - 5) ^ ^ mo \Ej\\ 



2{3-q) -p p{ \ \Ej\ niQ J 

2p(2-p)s^-{12p(g-l) + 8(3-g)}s-3p+12(3-g)-6(g-l)(3g-5) ^^ 

2{3-q)-p ' 

Consider the case when /^o = 0. From ( |112D we have = 0, and hence ( |129| ) implies that 
(5^ if > is equivalent to 

2p (2 - p) s^ - 4p (3g - 2) s + 3p - 6 (g - 1) (3g - 5) > , /^o = . (130) 

The solution of this inequality is 

A-p \2-p\ ^4-p |2-p| 

^ a — <^<^^ + — a — ' /^o = 0. (131) 

3 o 3 o 

Notice that for p = (|131|) coincides with the local stability condition proved in ref. ITO 



which we quoted in equation ( |103|) . Also notice that the only positive integers p and g which 
can satisfy ( |131| ) are p = 1 and g = 1, which is the case investigated in this paper. To 
have a sensible field theory p and g must be positive integers. Therefore the result - that 
stability alone restricts the choice among the class of theories defined by ( |105|) to just one 
case: p = q = 1 - must be considered as satisfactory. 

For each p and g obeying ( |131|) , the local stability for a sufficiently small /io 7^ follows 
by continuity of H4. The local stability condition 6^H > implies an upper bound on fiQ. 
For instance, for p = 1 and g = 1 we substitute equation ([75|) in ( p.29| ), which in the present 
notation reads 

H, = PmoN (^-\^] , (132) 

and obtain the local stability condition for the ground state 

2 oi+M , 3 1 + 14/3 ^ ^ 
^ -'TT^^+2TW^°' ^'''^ 
which must be satisfied for all real s. Solving for /3, we obtain 

/3<^(l + ^). (134) 

and then using (|77|) and (|79D 



cr<1.0, /io<A/2mo. (135) 

Thus a locally stable ground state for the theory investigated in this paper (p = g = 1) exists 
only if cr < 1. 
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Appendix 

In this appendix we develop a simple algorithm for the numerical determination of the 
functions needed for the evaluation of the functions and Vo-. As a first step we split the 
recurrence relation ( ^21) into three separate pieces to avoid multiple integration 

ry 

^nm{y,Cr)= dt (f)nmit,a) , (136) 

Jo 



(Pnm {y, cr) = -2 Vn-l,m {V, 0") 

[y + ncr + m + - f V 

n m 



+ y2J2 dt(pn-k,m-i{y -t,a) fki{t,a)] (137) 



k=0 1=1 



and 

fnm {y, cr) 



{y + na + m)^ — cr^ 




dt(pn-k,m-i {y - t, a) (t)k,i~i {t, a) . (138) 



It is easy to see that these equations are explicit recurrence relations, rather than integral 
equations. In order to convert these recurrence relations to a form digestible by computers, 
we divide the interval [0, y] into M pieces each of length x and use the trapezoidal rule to 
evaluate the integrals. As a result of the discretization we obtain a new set of functions 
^nm i^j^ cr)> 0nm i^j^ c^) and {xj, a) defined on the grid of points Xq = 0, Xi, . . . , xj^ = y, 
which converge to the true functions (pnm {xj, a), (pnm {xj, a) and fnm {xj, cr) as a; ^ 0. The 
corresponding recurrence relations are 

X / \ 

^nm i^j, ^) = Vnm i^j-l, ^) + ^ [Km i^j-l, (^) + 4>nm i^j, (^)) , (139) 



-T / \_ 1 /_ / ^ I { \ , ^ '^n,m-l[^jjCT) 

fnm K^ji^l — , s 2 ,1 Vn-l,m ^rj + \^Xj, (Jj + 



[xj + rid + m + - J 



n m j — 1 

+ ^J2J2Y1 ^n-k,m~i i^j - cr) Iki (a^i, (y) I (140) 

k=0 1=1 i=l 
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and 



fnm [Xj, = 7 ■ ■ ^2 ; V'n.m-l {^j, ^) + TT 0n,m-l i^j, 0^) + TT 7 

[Xj + na + m) — a-' \ ' I I a (cr + 1 j 

71 m J — 1 \ 

+ ^ mm ^n-Km-l (a^i - a^i, f^) (a^i, , (141) 

fc=0 i=l i=\ j 

where and j = 1, . . . jA/", n = 0, 1, ... , m = 0, 1, ... . The initial data for the 

recurrence process are 

^-i,m(a;i,f^) = 0, = 0, ^„^(0,a) = 0, (142) 

except for n = m = in which case 

^oo(0,^) = l, (143) 

and 

0„,_i(x,-,a) = 0, 0,^(O,a) = O, (144) 
except for n = 1, m = in which case 



The recurrence process starts with the evaluation of the cycle (|141|) — > ( |140| ) — > (|139|) for 



n = m = Q and j = 1. Then n, m and j are iterated until certain maximal values, jmax, 'T-max 



and mmax say, are reached, which are determined by the step-function in ([3T|). The result of 
the recurrence process is 7p^^ (xi, a) , . . . ,^„^ (a;j,„,,, cr), n = 0,..., n^ax, m = 0, . . . , mmax- 
In order to obtain an estimate of the difference between Tp^j^{xj, a) and Lpnm{xj,cr), the 
grid is refined by replacing A/" by 2Af and x by a;/2, and the recurrence process is re- 
peated. The resulting Tp^^{x2,<j) ,Tp^j^{x4,a) , . . . are compared with the previously cal- 
culated (xi, cr) (a;2, cr) , . . .. The process terminates when a certain specified ac- 
curacy is reached and the (xi, cr) ,(Pnm (^2, o") , . . . can be considered to be identical to 
<Pnm {xi, cr) , ipnm {^2, cr) , . . . . Notice that the numerical evaluation of the functions (pnm does 
not mean that we are solving the problem numerically. Rather, it means that the solution 
(|3l| ) is given in terms of nonstandard functions and that we have to teach our computer to 
obtain the values of these well-defined analytical functions. 
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b) The mass to bare mass ratios M/Nmp (solid line), Mqnn /Ntuq (dashed 



line) and M^qn /Islmp (dotted-dashed line) as functions of a for various hxed 



values of /io/mp. Dotted vertical lines correspond to the maximal values ol 
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minimum of the mass M/ JStuq = 2V2/3, which occurs at /i = 0.| 29 



2 a) The wave function ■0 (r) plotted as a function of the dimensionless variable 



7r for various values of the parameter a = /io/27, where 7 = (m^ — E"^ 



The function is scaled by an appropriate factor in order to become a dimen- 
sionless wave function of only two variables 7r and a. 'i'he latter depends only 



on the fundamental parameters of the model, such as the coupling constant g.\ 29 



b) The same as in figure^ for the potential function —(p (r) 29 
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3 a) The function R^j (r], C, s) as a function of s for a = 0.1 and r] = 2.448(5), 


C = 582.(3). The constants rj and C were determined from the boundary 


conditions (28) and (30) 
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4 'i'he constituent energies e/rriQ and e* /rriQ as functions of the coupling constant 
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(Jo for the maximum and the minimum of g respectively, according to the 


approximate formulae given in the text 



Figure 1: a) The mass to bare mass ratios Mn/Nttlq (solid line), Mqnn /Nitlq (dashed line) 
and Mmqm /NrriQ (dotted-dashed line) as functions of g'^N/Aiiml for various fixed values of 
Ho/rriQ. Dotted vertical lines correspond to the maximal values of g'^N/iTiml, the dotted 
horizontal line corresponds to the absolute minimum of the mass M/NniQ = 2V2/3, which 
occurs at /i = 0. 
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Figure 1: b) The mass to bare mass ratios M/Nuiq (solid line), Monn/Nuio (dashed hne) 
and MjqQjq /Ntrq (dotted-dashed line) as functions of a for various fixed values of fio/nio. 
Dotted vertical lines correspond to the maximal values of g'^N/ATcniQ {a = Ui), the dotted 
horizontal line corresponds to the absolute minimum of the mass M/NniQ = 2V2/3, which 
occurs at /i = 0. 



Figure 2: a) The wave function (r) plotted as a function of the dimensionless variable jr 
for various values of the parameter a = /io/27, where 7 = (ttt-q — E'^)^^'^ . The function is 
scaled by an appropriate factor in order to become a dimensionless wave function of only two 
variables 7r and a. The latter depends only on the fundamental parameters of the model, 
such as the coupling constant g. 



Figure 2: b) The same as in figure |^ for the potential function —(p (r). 



Figure 3: a) The function i?o-(?7,^,s) as a function of s for a = 0.1 and r] = 2.448(5), 
( = 582.(3). The constants rj and ( were determined from the boundary conditions ( pHj) and 

(H). 



Figure 3: b) The same as in figure ^ for the function Va {rj, s). 



Figure 4: The constituent energies e/rrio and 6*/mo as functions of the coupling constant 
g'^N/Anml at various fixed values of Ho/niQ. The dotted line indicates the value l/^/2 which 
is the boundary value deviding the two functions. 



Figure 5: g'^N/Airml as a function of a for various values of fio/rriQ (solid lines). The vertical 
dotted lines and the vertical dashed line give the values of ai and ctq for the maximum and 
the minimum of g respectively, according to the approximate formulae given in the text. 



29 



